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B — > (j)Ks CP asymmetries as a probe of supersymmetry* 
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Department of Physics, KAIST 
Daejeon 305-701, Korea 

We consider sparticle contributions to time dependent CP asymmetry in B — > <j)Ks- As for the 
gluino-squark loop, LR or RL insertion is more likely to give recently observed S^k < than LL 
or RR insertion. Neutral Higgs contribution does not change S^k very much due to the strong 
constraint from B s — > We also show correlations among related processes such as direct CP 

asymmetry in B — > X s ^ and B s -B s mixing, and discuss theoretical motivations for desired strength 
of flavor violation. 
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Time dependent CP asymmetry in B — > ^i^s- mode is 
written as 



4» 



_ r(u°(t) 



r(s (t) - 0^) - 

-C^k cos(Ara^) 



T(B*(t)^^K s ) 
- S^k sin(Am d t), 



where Am^ denotes the mass splitting between the two 
mass eigenstates composed of B° and 5 . Current data 
on the coefficients of the sine and cosine terms and their 
Standard Model (SM) predictions are summarized in Ta- 
ble |D It shows a 2.7 cr discrepancy between the measured 
value and SM prediction of S^k- This has become a hot 
issue in the high energy physics community recently. One 
of the reasons why physicists are focusing on5^ <pKs 
mode is that it is loop suppressed in the SM because it 
is a purely flavor changing neutral current (FCNC) pro- 
cess. This means that it is more sensitive to possible new 
physics effects which presumably arise at loop level, than 
those processes that have tree level contribution such as 
B — > J/i/jKs. In this work, we consider a possibility 
of filling this gap between the data and SM value, with 
effects from supersymmetry (SUSY). 

One major source of flavor violation in the general min- 
imal supersymmetric standard model (MSSM) is non- 
diagonal squark mass matrix. An off-diagonal element 
in the squark mass matrix in the super CKM basis can 
constitute a diagram leading to b — > sqq through gluino 





S(j)K 


C^K 


BaBarfl] 


+0.45 ± 0.43 ± 0.07 


-0.38 ±0.37 ±0.12 


Belle[2] 


-0.96 ± 0.501°;?? 


+0.15 ±0.29 ±0.07 


Average 


-0.15 ±0.33 


-0.05 ±0.24 


SM 


0.734 ± 0.054 





Avg. - SM 


-2.7a 


-0.2a 



TABLE I: Current status of time dependent CP asymmetry 
in B -> (j)K s . 
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FIG. 1: Gluino-squark loop penguin diagram. 



mediation, where q = u,d,s,c,b. An example of this type 
of QCD penguin diagram is shown in Fig. E] Here we 
adopt the mass insertion (MI) approximation which is 
one of the most convenient ways to analyzing this type 
of contribution. We consider four Mi's relevant to b — > s 
transitions, i.e., (^ L ) 23 , 0^)23, 0^)23, and (^ L ) 23 , 
one at a time. In addition to these four cases, we in- 
clude another one called RL dominance scenario |3(. In 
this case, CV 7 and Cg g are assumed to vanish at the 
scale, so the observed decay of B — > X s j must come from 
CV 7 , the chirality- flipped version of CV 7 . This unconven- 
tional scenario and the usual case can be distinguished by 
measuring the photon polarization. Therefore (6^)23 is 
forced to be a finite size to account for the data unlike the 
other four cases, and this corresponds to an extreme case 
with maximized SUSY effect in some sense. In the follow- 
ing, we turn on one of the four MI parameters at a time, 
and scan over it imposing constraints from B — » X s ^ and 
B s -B s mixing as follows: 

2.0 x 10" 4 < B(B X s7 ) <4.5 x 10" 4 , 
AM S > 14.9 ps _1 0. 

We impose a rather generous bound on B(B ^ X s j) to 
take into account theoretical uncertainties. The common 
squark mass fh and the gluino mass nig are chosen to 
be nig — m — 400 GeV. We use QCD factorization || 
in evaluating hadronic matrix elements for B — » (j)Ks- 
We do not assume any new physics in the B°-B mixing 
because sin 2(3 measurement from the B — > J/ipKs mode 
is well consistent with the SM fit. 
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(a) Allowed region for (5^ L )23 
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(b) S^x vs. C^x 
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(c) Sqk vs. AM, 
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(d) S<j> K vs. sin2/? s 



FIG. 2: (a) Allowed region on the complex plane of (^1^)23, and (b-d) correlations among observables for the LL insertion 
case. Different shades are used for different values of AM S . The black square corresponds to the SM case. The dotted line 
represents the upper bound on S^k at 1-a level. Lines in Fig. (a) are contours of S^k- 



First, let us consider the LL insertion case. Fig. [3 (a) is 
the region of (6^)23 consistent with B(B — > X s j) and 
AM S constraints. The solid lines are contours of S^k- 
The small square at the origin corresponds to the SM 
case, because there is no SUSY contribution there. The 
white space on the left is excluded by too small value of 
B(B X 8 i). The two holes around the Re(^ L ) 23 = 
line are excluded by too small AM S . In Fig. El (b), we 
plot predictions of Sqk and C^k from the allowed values 
of (Sf jL )23. The dotted horizontal line is the current 1-a 
upper limit on S^k- This figure shows that the minimal 
value of S(f,K in this case is about +0.5, and the change of 
S(f,K from the SM prediction is not significant. Although 
S^k is not very much affected, large effect in B s -B s can 
be expected. For example, we show the prediction of 
A M s in Fig. El (c). The SM value of AM S is about 16 
ps -1 , marked by the black square. We can read huge 
enhancement of AM S ~ 50 ps -1 is possible. Also, we 
show sin2/3 s in Fig. El (d), where 2/3 s is the phase of 



(B S \H^ B = 2 \B S ). The SM prediction of it is almost zero, 
but it can have any value between —1 and 1 here, which 
means that CP violation in B s -B s mixing can drastically 
change. 

The RR insertion case is almost identical to the LL 
insertion case. The only difference is the way B(B — > 
X s j) constrains (Sji R ) 23- The LL (RR) insertion mainly 
contributes to CV 7 (CV 7 ), and the dependence is such 
that 
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for LL, 
for RR. 



Because of this difference, (5^)23 plane does not show 
a region excluded by B(B — > X s j) that is present in 
Fig. 13 (a). However, this additional parameter space does 
not help shifting down S^k very much. Plots for this case 
are available in Ref. 

It is well known that the chirality-flipping LR and RL 
insertions receive enhancement by the factor of m^/m^, 
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FIG. 3: (a) Allowed region on the complex plane of (Sf jR )23, and (b-d) correlations among observables for the LR insertion 
case. The black square corresponds to the SM case. The dotted boxes represent current measurements at 1-a level. Lines in 
Fig. (a) are contours of S</>k- 



when they contribute to a (chromo) magnetic dipole op- 
erator. This is why they are strongly constrained by 
B — > X s j. For the same reason, they are more effective 
in modifying S^k and/or C^k than LL or RR insertion 
given the same magnitude of the MI parameter. We show 
the region of (Sf jR )23 consistent with B(B — > X s j) and 
AM S in Fig. 01 (a). As in the previous case, the solid 
lines are contours of Sqk, and the small square at the 
origin corresponds to the SM case. The annulus with 
radius ~ 10 -2 comes from the B(B — > X 3 j). Here the 
B(B — > X s j) constraint is so strong that AM S does not 
change very much from the SM value and play no role in 
constraining (5f jR )23- We show S^k and C^k from these 
values of (5f jR )23 in Fig. 01 (b). The dotted box is 1-a 
intervals of these observables, and we find some region 
where both S^k and C^k are within the box. This plot 
also shows a definite correlation between them. For S^k 
less (bigger) than the SM value, C^ K < (>) 0. Fig. 01(c) 
is the plot of B(B — > (j)Ks), with hatches on the region 



where B(B (j)K s ) > 1.6 x 1(T 5 . This region, excluded 
by B(B — > (j)Ks) constraint, corresponds to the hatched 
region in Fig. 01(a). Here we learn that significant por- 
tion of the (Sf jR )23 parameter space that is consistent 
with B(B -> X s i), is excluded by B(B -> <j>K s ). In 
Fig. 01 (d) is displayed the correlation between S^k and 
direct CP asymmetry in B — > X s ~f. For S^k less (big- 
ger) than the SM value, A b ^ sl > (<) 0. It appears that 
A^ S7 is outside the dotted box, but the experimental 
uncertainty is still large, so we cannot definitely conclude 
that this scenario is disfavored currently. 

We turn to the RL dominance scenario. Fig. 0] (a) is 
the region of (Sr L )23 constrained by B(B — > X s j) and 
AM S . Again, the solid lines are contours of S^k- In 
this case, we are assuming that SM contributions to CV 7 
and Cs g are somehow canceled by SUSY contributions. 
Because we are always involving SUSY contributions in 
CV 7 and C% g , there cannot be any point on the (Sr L )23 
plane that reduces to the SM case. The allowed annulus, 
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FIG. 4: (a) Allowed region on the complex plane of (£#l)23, and (b-c) correlations among observables for the RL dominance 
case. The black square corresponds to the SM case. The dotted boxes represent current measurements at 1-a level. Lines in 
Fig. (a) are contours of S</>k- 



centered at the origin, has radius ~ 10 -2 , which is fixed 
by B(B — > X s j). As in the previous case, the region 
consistent with B(B — > X s ^y) is always consistent with 
AM S as well. Predictions of S$k and C^k from the 
allowed region are shown in Fig. 21 (b). We can fit S^k 
to any value between —1 and 1. Moreover, big change in 
C^k from the SM value of zero, is expected in general, 
resulting in any value between —1 and 1. However, some 
region around the center is not covered. In Fig. 0] (c), 
we plot B(B — > 4>Ks), with hatches on the region with 
excessive B(B — > (j)Ks). Since there is only one weak 
phase in the b — > 57 amplitude, we have ^4 C ^ S7 = in this 
case. If we turned on another weak phase such as coming 
from (5^)23? we could have non- vanishing A^ S7 , but 
we do not pursue more detailed analysis here. 

We also analyzed the usual RL insertion case, keep- 
ing SM contributions to CV 7 and C% g . Because the SM 
prediction of B(B — > X s j) is already consistent with the 
data, there is less room for new physics in this case than 
in the RL dominance scenario. Nevertheless, we can get 
sizable change in S^k, satisfying B(B — » X s j). Plots for 
this case can be found in Ref. Q. 

These results of analysis depend on the sparticle mass 
scale. For a fixed value of MI parameter, SUSY effects 
get enhanced as the sparticle masses decrease. In Figs. 03 
we show possible values of S^k as functions of the gluino 
mass, rrig. In these plots, we fix the ratio x = m~/m 2 to 
be 1, and scan over one of the four MI parameters impos- 
ing constraints from B(B — » X s ^) and AM S . Figs. 0(a) 
and (b) are for the chirality-preserving insertions, LL 
and RR, respectively. Here we see a clear tendency of 
SUSY effects decoupling as the gluino mass increases. 
Note that our analysis shown in Figs.[2j0is for rrig = 400 
GeV. The way to getting S^k consistent with the data 
is to imagine the gluino mass ~ 250 GeV, which is close 



to the current experimental lower bound. On the other 
hand, a chirality- flipping insertion, LR or RL, shown in 
Figs. El(c) and (d), results in a range of S^k that is in- 
dependent of rrig. This is explained in the following way 
in the LR insertion case for example. The dependence 
of B — > X s j on the SUSY parameters is only through 
the ratio of {^LR)^/ m g^ which precisely is the param- 
eter that S<f>K also depends on. The set of (^LR)^/ m g 
allowed by B(B — > X s j) is independent of rrig, and so 
is S(f>K- The RL and RL dominance cases are explained 
likewise. In fact, this type of constant dependence would 
be true for LL or RR insertion case as well, if we al- 
lowed arbitrarily large size of (5 ( [ L )23 or ( 5^)23 and if 
we discarded AM S constraint. 

In the QCD factorization, divergences coming from 
endpoint of integral J dx/x are regulated in the form of 
X A(H ) = (1 + QA(H)e^ A ^)Mm B /A h ), where A h ~ 0.5 
GeV is the scale in hard-scattering. Discussion so far 
has been carried out without considering any uncertain- 
ties coming from this regularization, but we have to scan 
over the parameters Qa,h and (fA,H properly to take care 
of these theoretical uncertainties. In Fig. El (a), we show 
a curve of S^k as a function of the gluino mass. In this 
plot, we fix x = 1, and turn on only the RR insertion with 
the value (5$ 3 ) RR = 0.534 -0.856z, which minimizes S^k 
at rrig = 200 GeV. This curve is broadened into bands 
shown in Fig. El (b), if we begin to scan over the afore- 
mentioned parameters. In general, the uncertainty in the 
annihilation from qa dominates over that in hard specta- 
tor scattering from qh- The black band is for the range 
^ QA: Qh ^ 1 5 and this is what the authors of Ref. p 
recommend to use. On this band, the tendency of decou- 
pling SUSY effects with increasing rrig is still maintained. 
However, for the case of < Qa,Qh < 5, which is dis- 
played in light gray, estimated theoretical uncertainties 
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FIG. 5: Possible range of S^k as a function of the gluino mass. Different shades are used for different values of AM S . The 
dotted lines represent the current bounds on S^k at 1-a level. 



are so large that almost every value of S^k between —1 
and 1 is predicted, for any gluino mass between 200 GeV 
and 600 GeV. Ref. Q claims that they get S^k < from 
LL or RR insertion with m~ g = 350 GeV, disagreeing with 
us. We suspect this is due to the difference in choosing 
intervals for qa and qh- They use < qa < 8, and from 
Fig. [S| (b), it is obvious that S^k will nave every value 
between —1 and 1 for this interval. 

Let us think of theoretical motivations for (Sf jR )23 or 
($lr)23 ~ 10 -2 . In many alignment models using flavor 
symmetry, there remain small off-diagonal elements in 
the squark mass matrix in the super CKM basis. They 
can cause (S^ll or (S^rr ~ 10 -2 with arbitrary com- 
plex phase. This chirality preserving insertion can lead 
to an induced LR or RL insertion, 
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^rY^ - (^Ll)23 X 
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provided tan /3 is large enough so that \i tan /3 ~ 30 TeV. 
This kind of double insertion mechanism has been used in 
explaining both ex and c'/ck from a single CP violation 



source of (£f L )i2 8]. We also constructed a model that 
naturally gives the desired amount of LR or RL insertion 
using intersecting D5 branes p. 

Since QCD penguin contributes to 5 ^ sqq for 
q = u, 5, c, b at equal strength, one may worry 
about what happens to other decays such as B — > 
r] f Ks, ttK, K + K~ Kg. There are other 4-quark oper- 
ators which contribute only to these decays not affect- 
ing B — > <j)Ks, and changes in these decay modes are 
not definitely predictable in the present study. On the 
other hand, there is a nice mechanism which enables us 
to control B — > VP and B — > PP modes independently. 
Parity invariance tells us that the transition amplitudes 
of these modes depend on the new physics contribution 
in such a way that B — > VP (PP) mode depends on 
EJCf P + (-) Cf p ] #i, where Cf p and Cf p are Wilson 
coefficients with opposite chiralities, coming from new 
physics. Hence we can control B — > rj'Ks and B — > (j)Ks 
separately, and only the latter will change if we imagine 
a situation where (Sf jR )23 — (3rl)23i f° r instance. If this 
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FIG. 6: Predictions of S^k from (S^rr = 0.534 — 0.856z as functions of the gluino mass, (a) with Xa — Xh = 0, and (b) 
with varying qa^Qh- The dotted lines represent the current bounds on S^k at 1-a level. In Fig. (b), different shades are used 
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is the case, predictions of other observables such as C^k 
and A^ S7 will be different than the four cases we con- 
sidered here. After all, it requires further study to see 
changes in other decay modes in more general cases. 

From now, we consider the neutral Higgs mediated con- 
tribution to B — > (j)Ks- A typical diagram is shown 
in Fig. [7| This diagram contributes to B s — > /i+ \i~ 

i 



b 




s 



FIG. 7: Higgs penguin diagram. 

as well as to B — > (j)Ks if we replace the ss pair with 
Once we impose the upper bound from CDF Q, 
B(B S -> < 2.6 x 10~ 6 , we find that S^k > 0.71. 

This argument applies not only in minimal flavor viola- 
tion scenario, but also in general MSSM 6]. As a result, 
Higgs exchange does not cause substantial change in S^k- 
In conclusion, we analyzed SUSY contributions to time 
dependent CP asymmetry in B — > (j)Ks- As to the 
gluino-squark loop, negative S^k is more likely to come 
from LR or RL insertion than from LL or RR insertion 
if the sparticle masses are not close to the current exper- 
imental lower bound. However, we may have chances to 
look for LL or RR insertion in B s -B s mixing. Correla- 
tions among S<f>K, C^k, ^cp S1 •> AM S , sin2/? s may help 



us discriminate among different possibilities. Constraint 
from nonleptonic decays such as B — > <j)Ks is becom- 
ing as severe as that from B — > X s j. Higgs mediated 
FCNC cannot explain S^k < 0- Finally, let us direct the 
reader to Ref. 6] for more complete discussion including 
dilepton CP asymmetry in T(45) decay. 
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